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ABSTRACT 

In this paper we present Spitzer Infrared Spectrograph spectroscopy for 14 intermediate-mass young 
stehar objects. We use Spitzer spectroscopy to investigate the physical properties of these sources and 
their environments. Our sample can be divided into two types of objects: young isolated, embedded 
objects with spectra that are dominated by ice and silicate absorption bands, and more evolved objects 
that are dominated by extended emission from polycyclic aromatic hydrocarbons (PAHs) and pure H2 
rotational lines. We are able to constrain the illuminating FUV fields by classifying the PAH bands 
below 9 /im. For most of the sources we are able to detect several atomic fine structure lines. In 
particular, the [Ne II] line appearing in two regions could originate from unresolved photodissociation 
regions (PDRs) or J-shocks. We relate the identified spectral features to observations obtained from 
NIR through submillimeter imaging. The spatial extent of several H2 and PAH bands is matched with 
morphologies identified in previous Spitzer /IRAC observations. This also allows us to distinguish 
between the different H2 excitation mechanisms. In addition, we calculate the optical extinction from 
the silicate bands and use this to constrain the spectral energy distribution fit, allowing us to estimate 
the masses of these YSOs. 

Subject headings: ISM: jets and outfiows, ISM: lines and bands, ISM: photodissociation region (PDR), 
stars: formation, stars: pre-main sequence, stars: protostars 



1. INTRODUCTION 

Despite significant progress in our understanding of 
star formation, the earliest phases under which stars form 
are poorly characterized. This is particularly the case for 
stars of intermediate to higher mass (> 2 Mq), for which 
the evolutionary timescales are much shorter and which 
are smaller in number, compared to low-mass stars. 

Cloud cores that can form intermediate-mass stars 
have masses of ~ 10^ — 10^ Mq and temperatures of 
10 — 20 K and app ear in most cases as radio-quiet 
( [Arvidsson et al.|2010[ ) . Submillimeter studies of infrared 
dark clouds (IKDCs) also show high column densities 
(10^^ cm~^) capable of supporting intermed iate to high- 



mass star- formation (IVasyunina et al. 2009j). We have 
been studying a number of star forming regions, which 
are in an early evolutionary stage, matching these re- 
quirements, based on a large-scale unbiased sample of 
cold cloud cores or IRDCs (see Section [2]). Within these 
regions embedded mid-infrared point sources have been 
identified. A subsample of these sources are good candi- 
dates for intermediate mass young stellar objects (YSO). 
To study their physical conditions and chemical charac- 
teristics, we conducted infrared spectroscopy on these 
sources. 

One of the first infrared spectroscopic observation of 
a nearby star-forming region, performed by the Infrared 
Space Observator y (ISO), was a study of Orions IRc2 



van Dishoeck et a l. 1998; Wright et al.|2000D . Due to the 

large aperture of ISO's Short Wavelength Spectrograph 
(SWS), the measured spectrum is a superposition of mul- 
tiple physical components. The IRc2 spectrum contains 
hydrogen recombination lines, ionic fine structure lines 
and broad UV-pumped polycyclic aromatic hydrocarbon 



(PAH) emission bands originating from a photodissoci- 
ation region (PDR). Thermally excited, rotational, and 
rovibrational lines of H2 are also observed and can be 
used to infer the coupling between the molecular gas and 
dust grains. The PAHs and the smallest dust grains in 
this region are photoelectrically heated by a central hot 
object. These observations show solid state absorption 
features, such as silicates, water, methane and CO2 ice, 
originating from a quiescent extended ridge. In addition, 
shock excited molecular hydrog en was observed for th e 



adjacent Orion peaks I and II (Rosenthal et al. J 120001). 
A detailed study of the molecular gas-phase features of 
CO2, C2H2 and HCN appearing in absorption for IRc2 
and to som e extent in emiss i on for the Orion peaks can 
be found in |Boonman et al.] (2003). 

In this work, we use dcita o btained by the Infrared 
Spectrograph (IRS; | Houck et'a l. 2004) on-board the 
Spitzer Space Telescope ( [WernCT et al.| 2004). We find 
similar spectral features and use them to estimate the 
physical properties of the YSOs in the observed star- 
forming regions. From the 9.7 /im silicate features the 
optical extinction can be calculated and used f or later 
spectral energy distribution (SED) fitting (see Hen ning] 
(2010) for a general review on silicates). By classifying 
FAH bands below 9 /im it is also possible t o constrain 
some pr operties of the irradiating FUV field (Peeters et 
al. 2002] |Tielens|2008| . Class A PAH spectra are typical 
for H II and compact H II regions, refiection nebula, and 
YSOs exposed to a similar UV radiation field strength. 
The typical UV fiux intensities are 10^ — 10^ Gq. Re- 
ferring to jTielens^ ( ,2008 ) , class A spectra are excited by 
sources with '7eff >TDOOO K. This class of PAH spec- 
tra have an asymmetric feature at 6.2 /im, an emission 
complex peaking at 7.6 /im and a third PAH band at 
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8.6 /im. These PAH features arise mostly from the C— C 
mode. Strong UV radiation fields and gas densities of 
n ^ 10^ — 10^ cm~^ are needed to form photodissocia- 
tionregions (PDR) associated with molecular gas (Hol- 
|lenbach & Tielens 1997). The gas densities and the 
strength of the U V radiation field can be estimated from 
the ratios of certain forbidden atomic lines (e.g., [Fe II] , 
[Si II] ) and the excited H2 lines which we re created in 



the deep er layers of the illuminated PDR (Kaufman et 
|al.| 2006] ). It is not only in PDRs that emission lines 



01 ionized atoms can be observed. J-shocks interacting 
with the envelope of young protostars leave hot dense 
gas heated up to 10^ K. In the region behind the shock 
front, the molecules are dissociated and the ionized atoms 
can be observed by strong infrared emission lines (e.g., 
[Ne II] 12.8 /im, [Fe II] 26 /xm, [Si II] 34.8 /im, [S I] 25 /im) 
dHollenbach fc McKee 1989). The [Ne II] atomic fine 
structure line intensity depends on the shock velocity. 
Therefore strong [Ne II] lines can be only observed in 
high velocity J-shocks {vs > 60 — 80 kms~^). When the 
gas is cooling in the relaxation region further downstream 
from the shock front, H2 is formed and can be observed 
as pure rotational lines. In the non-dissociated C-shocks 
the temperatures never become high enough to dissociate 
molecular material. Spectral indications for C-sho cks are 



strong molecular lines of H2, CO, H2O, and OH (|Kauf- 
man & Neufeld 1996). Most spectroscopic studies of 
massive star-lornimg regions have used data from 750 's 
SWS and Long Wavelength Spectrograph (LWS). The 
IRS on board the Spitzer Space Telescope presents a dif- 
ferent instrumental approach. It has a lower spectral res- 
olution, but provides a better spatial resolution and an 
improved sensitivity compared to ISO. The long slit spec- 
troscopy in the low resolution mode allows spatially re- 
solved observations of selected lines down to a resolution 
of 1.8'^ Nevertheless multiple emission or absorption fea- 
tures from different spatial positions are still covered by 
the IRS beam. Therefore an interpretation of the spec- 
troscopic results can be given only when supplemented 
with multi- wavelength imaging observations. In general, 
these observations have a better spatial resolution than 
IRS in the same wavelength order and can therefore be 
used to trace certain spectral features. For example, the 
bands from the S pitzer^s Infrared Array Camera (IRAC; 
[Fazio et al.||2004|) and the Multiband Imaging Photome- 
ter (MIPS; ?) can be used to determine the spatial dis- 
tribution of excited hydrogen, PAHs, and warm dust. In 
particular the 4.5 /im band is ut ilized to identify the so - 
called green and fuzzy features ( Chambers et al.|[2009 ). 
These features predominantly originate from the (0-0) 
S(9) H2 band, mostly a ttributed to shock excitation (e.g. 



De Buizer fc Vacca 2010 ) . Another excitation mechanism 
for the line can be fiuorescence, depending on phys- 
ical parameters (e.g. temperature, H2 column den sities) 



and dust properties (iBlack fc van Dishoeckll 19871). Sub- 



millimeter observations (e.g. SCUBA, ATLASGAL) can 
also be used to trace cold dust. The current observations 
with Herschel will close the gap between the mid-infrared 
regime and the ground based submillimeter observations. 
For the cold 20 K) dust cores found, e.g. in the 



ISOPHOT 170 /im serendipity survey (ISOSS; [Bogun 
et al. 1996), the Planck spectrum peaks in the Herschel 



The next section describes the target selection. In Sec- 
tion [3j we describe the observations and the data reduc- 
tion processes. Section |4] details the continuum observa- 
tions and the spectral features. Section |5] discusses the 
different sources and their morphology within the con- 
text of previous observations. 

2. TARGET SELECTION 



Krause et al.| (|2004|) detail a number of suitable sources 
detected with ISOSS. Due to its high sensitivity, the sur- 
vey was able to find sources away from the galactic plane. 
We selected several regions within this sample based on 
their distance, luminosity and dust temperature (12 K 
< Td < 22 K). Later Spitzer imaging observations iden- 
tified point sources within som e of these re fflons (Birk-1 



2006^ ,2007| [Hennemann et al. |[2008, 20^ 



mann et al. 

We conducted a spectroscopic follow up of these objects. 
All our targets are resolved bright point sources at MIPS 
24 /im. Based on their mid-infrared colors, these sources 
were good candidates for intermediate-mass YSOs. The 
basic parameters of the selected ISOSS targets are listed 
in Table [1] (top part). In addition to the coordinates and 
the distances, the stellar mass estimates and the masses 
of the associated cold cores are given. In the last column 
the key references for the particular regions are listed. 

The ISOSS J18364-0221 East, ISOSS J20153+3453, 
ISOSS J04225+5150 and ISOSS J23053+5953 regions 
have extended PAH emission and warm dust compo- 
nents, which can be identified in the Spitzer /IRAC and 
MIPS 24 /im images, respectively. In the MIPS 24 /im 
band ISOSS J18364-0221 West appears as an isolated 
infrared source with weak dust components. 

To investigate additional targets similar to those cho- 
sen from the ISOSS sample, we selected 24 jam sources 
associated with c old cloud core candidates from I Cham- 1 
bers et al. (2009) that have IRS spectra available in the 
Spitzer-SSC archive. They are associated with IRDCs 
located in the galactic plane (called "IRDC sources" 
hereafter). They are clearly resolved point sources in 
the MIPS 24 /im and 70 /im bands (Table [ij bottom 
part). Except for GlO.70-0.13, 1.2 mm observations 
show compact cores with niasses ranging from 74 to 
301 Mq (Rath borne et al.| |2006|). For some of these 
sources H2O maser emission was observed indicating on- 
going star formation processes (Chamb ers et al. 2009]). 
Where ava ilable, the SCUBA legacy data ( pJi l^rancesco 



et al.|2008D for these regions reveal submillimeter clumps. 



In contrast to the other IRDC sources the GlO.70-0.13 
sources are not located within the dark filament of the 
IRDC but show extended PAH emissions such as, for ex- 
ample, the ISOSS J20153+3453 and ISOSS J04225+5150 
region. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Photometry 

We used the IRAC and MIPS obs ervations described 
in Birkmann et al. (2006 20071) and 'Hennemann et al. 



2009, 200^ to nieasure the integrated fluxes of the 
ISOSS sources. For the IRDC sources we used data from 
the Galactic Legacy Infrared Mid-Plane S urvey Extraor- 



FACS bands. 



dinaire (GLIMPSE; [Benjamin et al.|2003| a nd the MIPS 
Inner Galactic Plane Survey (MIPSGAL; Carey et al. 



2005 
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TABLE 1 

Sources 



Source 


R.A.(J2000) 


Decl.(J2000) 


Distance^ 


M ^ 

ivicore 




Key Ref. 




( hh ! inrQ ! ss ) 


(cici:iiirQ:ss) 


(kvc) 


(Mo) 


(Mr,) 






ISOSS 


















J04225+5150 East 


04:22:32.3 


+51:50:30 


5.5 






5.4-9.0 


(BPhD) 


J18364-0221 East 


18:36:36.0 


-02:21:59 


2.2 


75 




- 


(B06), 


(H09) 


J18364-0221 West 


18:36:29.6 


-02:21:59 


2.2 


280 




0.6-2.0 


(B06) 


J20153+3453 East 


20:15:21.2 


+34:53:46 


2.0 


87- 


149 




(H08) 




J20153+3453 West 


20:15:20.6 


+34:53:53 


2.0 


87- 


149 


4.0-8.5 


(H08) 




J23053+5953 East 


23:05:23.6 


+59:53:56 


3.5 


200 




4.1-6.6 


(B07) 




GOlO.70-0.13 I 


18:09:57.8 


-19:48:17.6 


3.7 






5.3-12.0 


(P09) 




GOlO.70-0.13 II 


18:09:55.9 


-19:48:17.6 


3.7 






5.6-6.1 


(P09) 




G019.27+0.07 


18:25:58.6 


-12:03:58.0 


2.4 


113 




2.1-4.9 


(R06), 


(C09) 


G025.04-0.20 


18:38:09.6 


-07:02:30.0 


3.4 


276 




5.4-9.6 


(R06), 


(C09) 


G028.04-0.46 


18:44:08.2 


-04:33:17.5 


3.2 


80 




1.3-7.6 


(R06), 


(C09) 


G034.43+0.24 


18:53:20.6 


+01:28:24.5 


3.7 


301 




3.6-6.4 


(R05), 


(R06), (C09) 


G053.25+0.04 I 


19:29:31.9 


+18:00:40.0 


1.9 






3.6-6.4 


(R06), 


(C09) 


G053.25+0.04 II 


19:29:31.7 


+18:00:33.2 


1.9 






2.3-6.1 


(R06), 


(C09) 



i 



References. — (BPhD) Birkmann I(l200 7|>; fB06) | Birkm ann et al.l(l2006); fB 07) Birkman n et al.| (|2007|>; 

^^^^ ^ fi annet al.TlBbU8t; (HU 9) Mennemann ^ al.| ( |2UUy| ; (PUy) [Peretto| 

Ubl; (HU^) pathborne e t al. (2UU6) 



C09) Chambers et al. (2009); (H08) Hennemann 
& Fuller (2UU9 ); (KU5j Kathborne et al. (2UU^ 



Note. — ^ The distances are obtained from the references in the last column. The kinematic distance 
for GOlO.70-0.13 is based on CO data from 'Dame et al.| ( |200l] >. 

^ The submillimeter core mass given in the key references for this object, if available. The targets in 
ISOSS J20153+3453 are separated by just 10^^ and the submillimeter peak is centered on the western source. 
^ The masses of the point sources were estimated by the SED fitting tool of Robitaille et al. (2007) based on 
the mid-infrared luminosities and the optical extinctions derived from the silicate bands if present. Because 
of the missing IRAC counterparts, no mass for the central source was modeled for ISOSS J 18364— 0221 East 
and ISOSS J20153+3453 East. 



The daophot package in IRAF was used to perform this 
aperture photometry. The resulting SEDs are shown in 
Figure [T] and [2] The apertures and annuh used for the 
sky background estimations for the Spitzer photometry 
are given in Table |3] Aperture corrections were not ap- 
plied to the IRAC mixes due to the inhomogeneous and 
asymmetric-extended background emission. The initial 
calibration of the IRAC and MIPS BCD data results 
in uncertainties of 2% for IRAC (Reach 2005), 4% for 
MIPS 24 mum fEn^elbracht et al.J l2007D and 10% for 
MIPS 70 /im (Gordon et al. 2(W|r Cosmetic corrections 
and astrometric enhancements we re perform ed using the 
MOPEX software package (Mako voz fc Mar le an 2005). 
The final images were combined using the IKAF frame- 
work. The resulting accuracy for the aperture photome- 
try is estimated to be 7% for the IRAC and 10% for the 
MIPS 24 /im. For the MIPS 70 /im band, we used point 
PSF photometry to decrease the uncertainties introduced 
by flux non-linearities (Gordon et al. 2007). The result- 
ing uncertainty for MIPS 70 /im is estimated to be 20%. 
For the ISOSS J18364-0221 East region the MIPS 24 /im 
peak was not detected as a point source in the IRAC im- 
age and as such photometry was only performed on the 
MIPS exposures. IRAC sources were not detected at the 
position of the eastern submillimeter peak in the ISOSS 
J20153+3453 region. A cluster of IRAC sources were de- 
tected at the position of the western 24 /im source. The 
photometry was performed on the central object of this 
cluster; this is also the position on which the IRS slit was 
centered. The eastern and western targets were resolved 
at 24 /im, but could not be detected separately in the 
MIPS 70 /im band due to the large beam size. All data 
points from the eastern and western sources are shown in 
the same plot (Figure [l|, and the data points longward 
of 70 /im are treated asbelonging to both sources. 



450 /im and 850 /im observati ons were obtained with 



SCUBA ( [Holland et al. 
Maxwell t elescop e (JUM " 



1999[ ) at the James Clerk 
I for mos t of the targets by 
(120061 120071 and [Hennemann et al. 
IS05SJ(M5^+5150 East was not in- 



Birkmann et al. 

(2008, 2oog'r 

eluded in these publications, but existing data from 2001 
were reduced for the 450 /im and 850 jam bands. For 
ISOSS J04225+5150 we derived submillimeter fluxes of 
F450 = (9.2 ± 2.8) Jy and Fgso = (1-5 ± 0.3) Jy from the 
SCUBA bands. A detailed descriptio n of the reduction 
and an alysis process can be found in [Hennemann et al.[ 
(2008). The absolute photometric accuracy is estimated 
toEe 30% at 450 /im and 20% at 850 /im. 

The low resolution spectra are plotted into the SEDs 
in Figures [l] and [2] The differences between spectra and 
the IRAC /MIPS photometry can be explained as follows: 
Our sources are compacted but extended. The available 
relative spectral response function is only calibrated for 
point sources (see next the section). The spectroscopic 
observations were conducted with a different aperture 
compared to the IRAC and MIPS photometries. This 
results in a significant flux difference. 

3.2. Infrared spectroscopy 

The spectroscopy was performed with the Infrared 
Spectrograph (IRS) on board the Spitzer Space Tele- 
scope. The observations of the targets selected from the 
ISOSS sample were performed during the campaign "IRS 
spectroscopy of extremely young massive proto-stars" 
(ID: 30919 + 40569). The observation dates, the pro- 
gram ID, the number of observational cycles and the 
integration times can be found in Table [2j To cover 
the maximum possible wavelength interval (5.2 jam — 
38.0 /im) all four available low resolution channels (SL2, 
5.2 - 7.7 /im; SLl, 7.4 - 14.5 /im; LL2, 14.0 - 21.3 /im 
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Fig. 1. — SEDs for the ISOSS targets were compiled using IRAC, 
MIPS, and SCUBA data. The used apertures for the Spitzer pho- 
tometry are stated in Table Is] The bar in the bottom right of 
each plot indicates an error oi 30% (absolute photometric errors 
are 7%, 10%-20%, and 20%-30%, for IRAC, MIPS, and SCUBA, 
respectively). The low resolution spectra are overplotted (see also 
Figures [3] and [4]). The fourth plot shows the results for ISOSS 
J201534'^453 East and West. Only for the western source are 
IRAC counterparts detected. The photometric data points (IRAC, 
MIPS 24 /im) for the western source are plotted as diamonds, and 
for the eastern source as circles (MIPS 24 /im only). Longward of 
70 /im both the sources are not separately resolved. The resulting 
data points are plotted as circles. 



and LLl, 19.5 — 38.0 /am) and both high resolution chan- 
nels (SH, 9.9 — 19.6 jam; LH, 18.7 — 37.2 /im) were used. 
ISOSS J18364— 0221 East was the only exception; here 
the SL2 channel was skipped due to the non-detection 
of a point source in the IRAC bands. The integration 
times for the low resolution modes were chosen to achieve 
a signal-to- noise ratio (S/N) of 100, sufficient to detect 
faint absorption features. In order to detect faint lines 
above the relatively bright mid-infrared continuum, we 
also aimed for a signal to noise ratio of at least 100 in 
the LH module at 24 /im. The integration times for the 
SH module were optimized in order to reach a spectral 
line sensitivity for H2 S{1) at A = 17.03 /im, which is 
two times better than that of the S{0) line. Overall we 
aimed for a line sensitivity better than 3 x 10~^^ Wm~^ 
and 5 x 10~^^ Wm~^ for the short and long wavelength 
orders of the high resolution spectra. All Astronomical 
Observation Templates (AOTs) were performed in the 
standard staring mode. The science targets were placed 
at 1/3 and 2/3 of the slit length to obtain spectra at two 
nodded positions. The slit overlays are shown for one 
nod position in Figures [9l p!3l p^H} [Tsl [191 and [25 
slit length covers several ten to hundredthousanc . 
as indicated in those plots. 
The IRDC sources G025.04-0.20, G028.04-0.46, 



The 
AU, 



G034.43+0.24 and G053.25+0.04 I+II were taken in the 
short wavelength orders only. The observations were car- 
ried out with a minimum of 10 cycles. The good sampling 
results in a high S/N. 



10^ 
10^ 
10^ 
10^ 

10^ 

10^ 
10^ 

10^ 
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10^ 

10^ 

10^ 
10^ 
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10^ 

10^ 

10^ 
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H ^ ^ 
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AT - 
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H ^ ^ 

/ — ^ 
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A 


A 
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Fig. 2. — SEDs for the IRDC sources were compiled using IRAC 
and MIPS data. The bar in the bottom left of each plot indicates 
the absolute photometric error of 10% for the MIPS 24 /im band. 

3.2.1. Low resolution data extraction 

The low resolution data products were obtained in 
the droopres state processed by the Spitzer pipeline. 
For further data reduction a pipeline based partly 
on the SMART-package together with the spectral 
extraction tools from the FEPS Spitzer science legacy 
te am were used. A detaile d des cription can b e found 
in Bouwman et al. (2006) and Swain et al. (2008). 
During the reduction process the data products were 
flatfield and straylight corrected. The correction of 
rogue and bad pixels was done with a median filter 
(based on the irsclean routine) and visual inspection. 
The spectral extraction was done with a fixed 6 pixel 
wide aperture for the short-low (SL) orders and 5 pixel 
wide aperture for the long-low (LL) orders, referring to 
spatial widths of 10. 8^' and 25. 5^' in the slit respectively. 
The local background was subtracted using the two nod 
positions along the slit from a median combined series 
of exposures. After the extraction the spectrum was 
convolved with the relative spectral response function 
(RSRF). The RSRF for the IRS instrument is only 
defined for point sources. Due to the compact, but 
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TABLE 2 
Observation Log 



Source 


AOR-ID 


Obs. Date 




tramp (s), (cycles) 




Comments on 






(yyyy-mm-dd) 


SL 


LL 


SH 


LH 


Data Quality 


ISOSS 
















J04225+5150 East 


19235328 


2007 Mar 15 


6,(2) 


6,(2) 


30,(2) 


14,(2) 




J18364-0221 East 


23071744 


2008 Apr 26 


60,(3) 


30,(2) 


120,(6) 


60,(8) 


SL2 not observed, bright background emission. 








SL orders not extracted 


J18364-0221 West 


23072000 


2008 Apr 26 


6,(2) 


6,(2) 


120,(2) 


60,(2) 




J20153+3453 East 


19236864 


2006 Jul 04 


14,(4) 


14,(2) 


120,(3) 


60,(4) 


The western source dominates SL orders 


J20153+3453 West 


19237376 


2006 Jul 4 


6,(2) 


6,(2) 


30,(5) 


14,(5) 




J23053+5953 East 


19237888 


2006 Jul 31 


60,(3) 


6,(2) 


120,(2) 


60,(2) 


Flux leakage in SL orders 








Extended PAH band^ 



GOlO.70-0.13 I 
GOlO.70-0.13 II 
G019.27+0.07 
G025.04-0.20 
G028.04-0.46 
G034.43+0.24 
G053.25+0.04 I 
G053.25+0.04 II 



16943872 
16943872 
12091648 
17823744 
17823488 
17824512 
17825280 
17825536 



2006 May 1 

2006 May 01 
2005 Apr 18 

2007 May 08 
2007 May 02 
2007 May 02 
2007 Jun 14 
2007 Jun 14 



6,(1) 
6,(1) 
60,(2) 
60,(10) 
60,(10) 
14,(34) 
60,(10) 
60,(10) 



6,(1) 
6,(1) 
14,(2) 



6,(1) 
6,(1) 



14,(1) 
14,(1) 



very non-uniform background^ 

extended PAH bands^ 
extended PAH band^ 
extended PAH band^ 
extended PAH bands^ 



Note. — ^ The 11.2 fim PAH band is spatially extended over the slit. A polynomial fit is used to estimate the background 
contribution in this wavelength range. 

^ Because of the very non-uniform background the whole spectrum was extracted with a polynomial background fit. 

^ The 11.2 fim and 12.3 /im PAH bands are spatially extended over the slit. A polynomial fit is used to estimate the background. A 

polynomial fit is used to estimate the background contribution in this wavelength range. 



TABLE 3 

Photometry Apertures for IRAC and MIPS 



Source 



IRAC 



MIPS 





^Ap 


Rsky 


^Ap 


^Sky 


ISOSS 










J04225+5150 East 


2^' 






13'' 


J18365-022 East 








13" 


J18365-022 West 


3.2^^ 


h.2" 


13.2'' 


30" 


J20153+3435 ^ 


2" 






13" 


J23053+5953 


2" 


Q" 


Q" 


13" 


GlO.7-0.13 I 


4(^8 


10(^8 


10'.'4 


31'.' 2 


GlO.7-0.13 II 


4(^8 


10(^8 


10'.'4 


31" 2 


G19. 27+0.07 


4(^8 


12(^0 


10'.'4 


31" 2 


G25.04-0.20 


r!2 


14(^4 


10'.'4 


31" 2 


G28.04-0.46 


7'!2 


14(^4 


10'.'4 


31"2 


G34.43+0.24 


7'!2 


14(^4 


10'.'4 


31"2 


G53.25+0.04 I 




19^6 


10^4 


31"2 


G53.25+0.04 II 




19^6 


10^4 


31'.'2 


Note. — ^ Only the western source is detected 



within the IRAC bands. For MIPS at 24 /xm the 
same aperture size was used for the eastern and 
western sources. 

extended nature of the sources a proper absolute flux 
calibration cannot be achieved for the one-dimensional 
spectra. Since there is a flux offset between the SL and 
LL orders and the corresponding slits have different 
orientations, the one-dimensional spectra are shown in 
two separate plots in Figure |3] Furthermore, the fluxes 
are given in units of Jy with respect to the used aperture 
sizes. 



3.2.2. High resolution data extraction 

The high resolution data (SH: short-high orders; LH: 
l ong-high orders) were re duced with the c2d pipeline 
(Lahuis et al. 2006 2010[ ) starting with the rsc data 
products. Corrections tor bad and rogue pixels were 
performed, as was defringing for all spectral orders. 



The optimal spectral extractions were obtained using 
analytically derived cross-dispersion proflles that flt the 
source proflle and extended emission components. The 
background was determined from additional off-source 
AOTs for every source. We used the high-resolution 
data to obtain line fluxes and flux ratios. 



3.2.3. Limitations on the extraction process 

To analyze the effect of extended emission lines and an 
inhomogeneous background on the low-resolution spec- 
tra, we investigated the dependence of the sources' cross- 
dispersion proflles as a function of wavelength. Several 
spectra showed extended spatial proflles, in particular 
for the 11.3 /im and 12.8 /im PAH bands. These ex- 
tended PAH background contributions could lead to a 
background over-subtraction in the resulting low resolu- 
tion spectra. We compared different extraction methods 
for the low resolution spectra. We used a polynomial 
background flt for the FEPS extraction in the affected 
wavelength regimes. The second method was the PSF- 
extraction from the c2d pipeline, which is slightly un- 
dersampled for the low resolution case. Both extraction 
methods do not signiflcantly differ in the resulting over- 
all continuum, correcting for extended background line 
emissions. For some sources the c2d extraction returns 
a slightly better S/N for the continuum. In some cases 
the c2d pipeline did not detect some of the faint H2 rota- 
tional lines. In general the performance of fltting the SL 
order is better with the FEPS extraction. Based on these 
facts we chose to process sources which were affected by 
extended PAH emissions with the nominal FEPS extrac- 
tion and substitute the affected wavelength intervals with 
the polynomial background flt (see Table [2]). 

ISOSS J18364-0221 East shows bright and non- 
uniform background emissions at short wavelengths 
(5.2— 14. 5'^). No source proflle could be flt for these wave- 
length orders and as such no SL spectrum is presented 
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10 11 12 13 1415 
X [|im] 



30 35 



Fig. 3. — Low-resolution spectra for the ISOSS sources. The SL and LL orders are extracted with different aperture size. Therefore the 
long wavelength orders are more affected by contributions of extended emissions and show a clear offset compared to the short wavelengths. 
The low resolution, long wavelength orders for ISOSS J18365-0221 East and ISOSS J20153+3453 East are shown in Figure [i] Emission 
lines are indicated by lines, while absorption bands are underlined in light gray. 



in Figure [S] For the ISOSS J20153+3453 region two tar- 
gets were observed centered on the two MIPS peaks at 
24 jam. The separation of both MIPS peaks at 24 /im 
is 10'^ The sUt positions of the SL orders overlap for 
the AOTs taken on both these targets. For the eastern 
target in this region we did not detect a point source in 
the IRAC bands, but instead detect a diffuse bulge of 
emission. The independent source profile could not be 
obtained for the SL orders of the eastern target since the 
slit is dominated by the western source. The short or- 
ders were therefore not extracted for the one-dimensional 
spectrum. The long wavelength orders (LL) are not af- 
fected as these slit positions do not overlap. 

The high-resolution spectrum could also not be ex- 
tracted for this target in the SH orders as the cross 
dispersion profile could not be fitted with the absence 
of a point source in SH orders. For J23053+5953 East 
the slit was centered on the peak of the submillimeter 
clump. The IRAC point source at this position was only 
partly covered by the slit, which results in a significant 
flux loss. We can therefore only qualitatively discuss the 
spectrum below 14.5 /im. G019.27+0.07 is affected by 
a very non-uniform background with bright background 
emissions. The whole spectrum was extracted with third- 
and fourth-order polynomial background extractions for 
the SL and LL orders, respectively. It is shown in Fig- 
ure [7| together with the averaged background spectrum 
at the nominal nod position. 



3.2.4. Spatial line analysis 

We also investigated the spatial extent of single lines 
over the low resolution slits (SL,LL). A one-dimensional 
spectrum was extracted for every spatial pixel position 
over the slit width. This extraction was performed with- 
out a background subtraction. The spectra for the two 
nod positions were extracted separately and were not me- 
dian combined. The line flux was calculated from a fitted 
Gaussian profile after the underlying continuum contri- 
bution had been removed using a second-order polyno- 



mial. These results are presented in Section [46 
4. RESULTS 

4.1. Continuum observations and spectral energy 
distribution 

The SEDs were compiled from IRAC, MIPS and 
SCUBA bands and are presented in Figure [l] 

Using these data we calculated the bolometric lumi- 
nosities using: 

L = And^ J dvF^ 

This assumes a spherically symmetric distribution of the 
luminosity. Of course, other non-spherical configurations 
are possible; for example, the majority of the observed 
luminosity escapes via outflow cavities. For this reason, 
the values presented in Table |4] are affected by large un- 
certainties. 
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TABLE 4 

Physical parameters and Spectral features 



Source 


Distance 


Lboi 


PAH 


Ices S 


z Molecular absorptions 


Silicate 


TT h 




Atomic lines^ 




























IRc2(A) 


0.45 




/ 


H2O, 


CO2, 


C2H2 


/ 


/ 


Ne II 


1 , 


[Ne III] , [P III] , 




















Fe II 


, [S I] 


[S III] , [Si II] 


ISOSS 


























J04225+5150 East 


5.5(B) 


5000 


/ 


CO2 






/ 


S(0)-S(4), S(5) 


Ne II 


] - 


[Fe II] , [Si II] 


J18364-0221 East 


2.2 


70^ 


/ 


_ 






_ 


S(0) 


[SI] , 


[S 


n] , 


[Fe II] 


J18364-0221 West 


2.2 


550 


- 


CO2, 


H2O, 


CH4 


/ 








J20153+3453 East 


2.0 


920 


/ 


CO2 






/ 


S(l) 


Fe II 




Si II 




J20153+3453 West 


2.0 


920 


/ 








emiss. 


S(l), S(2) 


Fe II 




Si II 




J23053+5953 East 


3.5 


2300 


/ 


CO2 






/ 


S(0)-S(7) 


Fe II 




Si II 


, [S I] 


GOlO.70-0.13 I 


3.7 


39^ 


/ 


H2O, 


CO2 




/ 


S(l) 


Ne II 


] - 


[Si II 


] 


GOlO.70-0.13 II 


3.7 


18^ 


/ 








/ 


S(l) 


Ne II 


] , 


[Si II 


1 


G019.27+0.07 


2.4(C) 


9^ 




C02^ 








S(0)f, S(lY no hi- 


res 




G025.04-0.20 


3.4(C) 


4ie 




H20, 


NH+ 


CH4 


/ 


no hi- 


res 




G028.04-0.46 


3.2(C) 


39^ 


/ 


H20, 


NH+ 


CH4 


/ 


S(l) no hi- 


res 




G034.43+0.24 


3.7(D) 


25^ 


/ 


H20, 


NH+ 




/ 


S(2), S(4)-S(7) no hi- 


res 




G053.25+0.04 I 


1.9(C) 


5^ 


/ 


H20, 


NH+ 




/ 


S(2)-S(7) no hi- 


res 




G053.25+0.04 II 


1.9(C) 


4e 


/ 


H20, 


nhJ 




/ 


S(4), S(5) no hi- 


res 





References. — (A) |van Dishoeck et aL] ( |1998t ; (B) [Birkmann | ( |2007t ; (C) [Rathborne et al.| ( |20Q6t 



Note. — 

^ The bolometric luminosities were estimated using the SEDs shown in Figure [l] 

^ The presence of pure rotational lines in the high or low resolution spectra are listed in this column. 

^ Atomic lines which could not be observed due to missing high resolution spectra, were labeled as not available (no hi-res). 
^ Due to the missing IRAC counterparts the calculated luminosity for ISOSS J 18364— 0221 East represents just an lower limit. 
^ The luminosities for the IRDC sources are calculated with the IRAC and MIPS bands only. 
^ Only observed in the background spectrum. 

The low resolution spectra extracted from the ISOSS- 
sources are presented in Figure [Sj The spectra from 
the IRDC sources and the long wavelength orders from 
ISOSS J23053+5953 East and ISOSS J18364-0221 East 
are shown in Figure |5] IVIost spectra in our sample are 
dominated by broad emission bands at shorter wave- 
lengths. Beyond 18 /im the continuum is well defined 
and shows a steeper slope compared to the shorter wave- 
lengths. For ISOSS J23053+5953 East the IRAC obser- 
vations reveal a cluster of objects that are not resolved 
at the longer wavelengths. The contribution from the 
clustered sources can be seen in the higher fluxes below 
9/im. Clustering can also be seen in ISOSS J04225+5150 
and ISOSS J20153+3453 West. 

4.2. Hydrogen emission 

The pure rotational lines {v' — v^') = (0 — 0) detected in 
our sample range from S{7) to 5'(0), with corresponding 
excit ation energies of Es^ I k = 510 K to Es^/k = 7197 K 
(e.g., [ Rosenthal et al.||2000D . 

For 12 out of 14 sources we observed rotational II2 
lines. Besides thermal excitation, other mechanisms such 
as shock excitation or PDR contributions at the atomic- 
to-molecular interface must be considered (see Section 
|5|. Except for three IRDC sources, the low J < 2 lines 
with low excitation energies < 1015 K which can trace 
warm gas were detected. Only for one source were all 
the (0-0) transitions present, from S'(O) to 5'(7). Due to 
the low S/N ratio in the low resolution spectra the S{0) 
line can only be identified in the high resolution spectra. 
The observed rotational lines are listed in Table ffl The 
high resolution spectra reveal the presence of the 5(2) 
(1-1) line for ah ISOSS sources apart from the ISOSS 
J20153+3453 targets. 
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Fig. 4. — Low-resolution spectra for ISOSS J18365 East and 
J20153+3453 East are shown. Only the long wavelength orders 
(LL) were extracted for these sources and plotted. 

For ah targets with detected S{2), S{1) and S{0) lines, 
the line ratios S{2)/S{1), S{2)/S{0) and ^(l)/^(0) are 
below 1. 

4.3. PAH features 
The PAH bands were identified with the features tab- 



ulated in Draine & Li L2007). The 6.2 /im and 7.7 /im 
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Fig. 5. — Low-resolution spectra for the IRDC sources. The spectra for G10.7— 0.13 and G19.27+0.07 are shown in Figures [6] and [T] 
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Fig. 6. — Low-resolution spectra for GIG. 7— 0.13 I and II are shown. 
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features originate from pure CC stretching and CH in- 
plane bending. They, alon^ wit h the 8.6 /im CH in-plane 
vibration mode fTielens"2008'), can be clearly detected 
in ISOSS J04225+5150, ISOSS J20153+3453 West, and 
GOlO.70-0.13 II (see Figure[24|. These bands are also ob- 
served in the background spectrum of G019. 27+0.07, but 
not in the source spectrum. For these targets the first 
feature is centered at 6.21-6.22 /im and the PAH emis- 
sion complex around 7.7 //m is asymmetrical with a peak 
between 7.57 and 7.64 /im. In each case, a secondary, 
weaker component is embedded in the red flank of this 
profile. Based on the profile shape and peak position, we 
conclude that class A PAHs are dominant in these regions 
dPeeters et al]|2002[ ). The 11.2 jam and 12.8 /im features 
can also be clearly identified in the spectra from these 
sources. These bands origin ate from CH out-of-plane 



bending modes ( Tielens|2008 ). In addition, several PAH 
bands between 16.2 and 18.2 /im are present, although 



these features blend into a plateau (see Figure 



panel). Similar plateaus have been obse rved with ISO for 



24 



right 



some YSO s and compact H II regions ( [Van Kerckhoven 



|et al.|200"Q ^ . Another broad PAH band is clearly observed 
between 18.9 and 19.0 /im for ISOSS J04225+5150 East, 
ISOSS J20153+3453 West and GOlO.70-0.13 II. Tielens 
( [2OO8J argued that this feature could arise from highly 
ionized PAHs. 

The 11.2 /im band can appear as an extended bright 
background emission which is not associated with the 
source itself. This is the case for ISOSS J18364-0221, 
ISOSS J23053+5953, G028.04-0.46, G034.43+0.24, and 
G053.25+0.04 I+II. The same effect is also present 
for the 12.8 /im PAH band for G028.04-0.46 and 
G053.25+0.04 I. No PAH bands were detected in the 
background subtracted spectrum from GOlO.70-0.13 I, al- 
though the 6.2 /im feature was observed as an extended 
background emission. No PAH features were found in 
the spectrum from ISOSS J18364-0221 West. 

4.4. Silicates and molecular absorptions 

We were able to identify several absorption bands in 
our sample. The 9.7 /im silicate absorption feature 
was found in most sources of our sample. For ISOSS 
J20153+3453 West a broad emission component appears 
between 8 and 13 /im instead of the silicate absorption. 
We identified absorption bands by water ice (6.0 /im), 
NHJ (6.85 /im), methane ice (7.67 /im) and carbon diox- 
ide ice (15.2 /im). Those absorption bands and their op- 
tical depths are listed in Table [5] The IRDC sources, 
only taken in the short wavelength orders, are not acces- 
sible for CO2 ice detections. To calculate these optical 
depths, we fitted the spectral continuum by using mul- 
tiple blackbody components and removing possible PAH 
contamination. This was accomp l ished with the PAHFIT 
routine described in Smith et al. (2007). We calculated 
the op tical depth r for the ice features as in |Quanz et al. 
([2QQ7|) using 



In 



06s 



Here, Fmin{^) is the observed flux of the absorption fea- 
ture and FcontW is the fitted continuum flux. In addi- 
tion to the removal of PAH contributions we utilized the 
PAHFIT routine to calculate rg.y for the silicate absorp- 
tion. The routine is fitting the underlying continuum 



via multipl e blackbody fits and using the silicate profile 



from Kemper et al. 



see jSmith et al. 
culated from rg 



(For a detailed description 
The optical extinction Ay cal- 



is tabulated in Table |5| 



4.5. Ionic fine structure lines 

Several fine structure lines were detected in the high 
resolution spectra. The identified lines are listed in Ta- 
ble [4] and plotted in Figure |8] after application of a sec- 
ond order baseline fit. The iron line [Fe II] can be de- 
tected at 25.99 jam or 35.35 jam. The [Si II] appears at 
34.82 /im. The sulfur line [S I] (25.25 /im) is only de- 
tected for two sources in our sample. The neon emission 
[Ne II] (12.81 /im) always appears on top of a broad PAH 
feature. 

4.6. Spatial extent of lines 

The spatial line fits for the ISOSS J20153+3453, ISOSS 
J04225+5150 and the ISOSS J23053+5953 regions are 
shown in Figures [TT] , [15) and [22) respectively. In these 
figures the relative line fluxes are shown as a function 
of their spatial distributions. The positions for the line 
extraction are overlaid for the short (SL) and long (LL) 
wavelength orders and the two different nod observations 
(in red and blue). For the LL slit positions we indicate 
only the overlap of the two nodded observations on the 
IRAC map and not the whole slit length. On top we 
present the spatial line fits for the long wavelength or- 
ders (14-37 /im) for the S{1) 0-0 H2 line and the 16.45 /im 
PAH feature. The lines were fitted with a second-order 
polynomial for local background estimation and a Gaus- 
sian for the line fit. The bottom plots present the line 
fits for short wavelength orders (5.5 — 14.5 /im) for PAHs 
and molecular hydrogen lines S{b) and S{3). 

For the ISOSS J20153+3453 region we find similar spa- 
tial profiles of the 6.2 and 11.2 /im PAH features, which 
both peak slightly eastward (~ 4 — 5'^) of the western 
MIPS maximum at 24 /im. The line fluxes from both 
these PAH features decrease more quickly toward the 
east than toward the west. The relative line flux distribu- 
tion of these features is similar towards the western end of 
the slit, but, unlike the 11.2 jam feature, the 6.2 jam PAH 
band cannot be detected eastwards of the eastern MIPS 
source at 24 /im. The two molecular hydrogen lines 5'(3) 
and S{5) are not detectable at the western 24 /im MIPS 
peak location, but instead from two extended features to 
the east and west of this position. 

For ISOSS J23053+5953 East we detect an extended 
PAH feature at 11.2/im with almost constant line flux 
over the entire slit. Although the 6.2 /im PAH feature 
is not detected in the whole slit, it does have a line flux 
distribution similar to the 11.2 /im feature in the western 
end. The S(5) H2 feature is detected to both the east and 
west of the MIPS peak at 24 jam. The S(3) H2 line is also 
detected at the MIPS peak position. The 16.45 jam PAH 
feature is only detected in the vicinity of the MIPS 24 /im 
source , whereas the S{1) hydrogen line is extended over 
more than 180^' in this region. 

The spatial line fits for the ISOSS J04225+5150 region 
are shown in Figure [15) The line fluxes for the PAH fea- 
tures at 6.2 /im and 11.2 /im and the 5'(3) H2 line peak at 
the position of the brightest source observed in the IRAC 
and MIPS bands. The decrease of the line fluxes matches 



10 



Pitann et al. 




9 10 11 12 13 1415 20 25 30 35 

X \\im] 

Fig. 7. — Spectra for G19.27-0.07 are shown. The source spectrum is extracted with a polynomial background fit and is plotted in black. 
The source background spectrum is plotted in gray. 



TABLE 5 
Optical depths. 



Source 




^NHt 




^C02 


'^Silicate 


A ^ 


ISOSS J04225+5150 East 








0.342 ±0.012 


0.57 ±0.03 


7.9 


ISOSS J18364-0221 West 


0.256 ±0.017 


0.153 ±0.010 


0.055 ±0.009 


0.51 ±0.04 


1.38 ±0.07 


19.2 


ISOSS J23053+5953 East 








0.34 ±0.11 


4.8 ±0.3 


67 


GlO.70-0.13 I 


0.34 ±0.05 


0.24 ±0.04 




0.51 ±0.06 


1.60 ±0.08 


22 


GlO.70-0.13 II 










1.63 ±0.09 


23 


G025.04-0.02 


0.49 ± 0.05 


0.399 ±0.017 


0.118 ±0.013 




2.5 ±0.2 


35 


G28.04-0.46 


0.84 ±0.08 


0.33 ±0.04 






1.01 ±0.06 


14.1 


G34.43+0.24 


1.14 ±0.14 


0.44 ±0.04 


0.28 ±0.05 




1.7 ±0.09 


24 


G053.25+0.04 I 


0.348 ± 0.027 


0.174 ±0.012 


0.081 ±0.010 




0.267 ±0.013 


3.7 


G053.25+0.04 II 


0.43 ± 0.08 


0.80 ±0.08 


0.18 ±0.04 




0.47 ±0.03 


6.6 



Note. — ^ To calcula te the si licate optical depth T9.7 after the correction for possible PAH contamination we 
used the PAHFIT routine (Smithet al. 2007 ). ^ The op tical extinction A^ was calculated from the silicate bands 
using the extinction model trom [Hobitaille et al.| ( |2007| >. 




Fig. 8. — Fine structure lines as observed in the high-resolution spectra. 
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TABLE 6 
Line flux ratios for [Ne II] 



Source 


-^[Ne II] 


-^[Ne II] 


-^[Ne II] 


-^[Ne II] 


-^[Ne II] 


^S(2) 




^S(O) 


-^[Fe II] 


-^[Si II] 


ISOSS 












J04225+5150 


1.04 


0.63 


0.64 


0.22 


0.04 


GOlO.70-0.13 I 


2.8 


0.36 


0.04 




0.10 


GOlO.70-0.13 II 


2.3 


0.85 


0.44 


0.25 


0.5 



the morphology observed in the IRAC data. The slit 
for the LL observations partially covered the elongated 
submillimeter structure observed by SCUBA. While the 
16.45 /im PAH and the H2 S{1) line fluxes decrease in the 
near vicinity of the central source, the H2 line appears 
to be almost constant through the western submillimeter 
peak. 

5. DISCUSSION OF INDIVIDUAL SOURCES 

All the sources exhibit a complex morphology in the in- 
frared. Therefore, their spectra must be reviewed in the 
context of previous observations. The different sources 
will be discussed individually, sorted by their evolution- 
ary stage. The first sources are the most evolved ones, 
with indications for PDRs, followed by young, deeply 
embedded sources with several absorption features. 

5.1. The ISOSS J20 153^ 3453 region 




Fig. 9. Overlay of the IRS slits for ISOSS J20153+3453. The 
IRAC 8 /im (left) and MIPS 24 /im (right) images are shown in 
false colors. North is up and east is left. The slit positions for the 
IRS spectroscopy are overplotted for the low-resolution spectra in 
green (short wavelengths, SL) and red (long wavelengths, LL). The 
white overplots represent the high-resolution AOTs. The first nod 
position is indicated in the left panel and the second is in the right 
panel. The source is indicated by the dark magenta box. 

The ISOSS J20153+3453 regiom located 2.0kpc away, 
was studied by Hennemann et al.| ( |2QQ8J . In the MIPS 
observations at 24 /im, the clump is resolved into a 
eastern and western components with a separation of 
10". A single submillimeter clump was detected close 
to the eastern source by SCUBA. A dust temperature of 
Td = 15.3 - 17.0 K and a gas mass of M = 87 - 149 M© 



RA(2000) 

Fig. 10. IRAC composite image of the ISOSS J20153+3453 
region: the 3.6 /im channel is colorized in blue, 4.5 fim in green, 
and 8.0 /im in red. The contours of the SCUBA observations at 
450 /im are shown in yellow. The two MIPS point sources at 24 /im 
are indicated by the red box points. They are used for the IRS 
pointing. At the position of the eastern MIPS sources and the 
SCUBA source an elongated "green and fuzzy" feature can be seen. 

were calculated for the clump. Based on the IRAC col- 
ors, the MIPS source at 24 /im associated with the sub- 
millimeter peak was assumed to be a class O/I object. 
The mass reservoir in this region contains enough gas 
to form a massive star, but due to the multiplicity of 
NIR and IRAC sources the interpretation of the further 
evolution in this region is not straightforward. Another 
MIPS 24 /im source, surrounded by extended emissions, 
is located on the northwestern limb of the submillimeter 
clump. 

The IRS spectrum of the western MIPS 24 /im source 
is dominated by several PAH bands. The broad emission 
feature between 8 and 13 /im is typical for small, tran- 
siently heated dust grains distributed over a wide spatial 
range around the central sources. Photoelectric heating 
of the small dust and PAH grains accounts for thermal 
coupling of the grains to the gas, hence a certain amount 
of the observed molecular hydrogen must be thermally 
excited. The smooth spatial distribution of the 6'(T) H2 
line represents such thermally excited gas (Figure 11). 

As for the ISOSS J23053+5953 region, the detected 
[Si II] and [Fe II] lines and the spatial distribution of the 
S{5) and S{3) H2 emission can be attributed to J-shock 
excitation. The elongated "green and fuzzy" structures 
on the IRAC composite images fit this picture as well (see 
Figure[lO|. In contrast to ISOSS J23053+5953, the PAH 
features m the spectrum of ISOSS J20153+3453 West are 
identified as class A according to ^Peeters et al. (2002 ) (see 
Section fl|. Therefore, UV field strengths are sufficient to 
form a FDR. The 18.9 zim fe ature could represent highly 
ionized cationic PAHs ( [Tiel ens 2008 ) and give further ev- 
idence for a strong UV radiation. The detected [Si II] and 
[Fe II] lines are found not only for J- sh ocks; they are 
known as minor c oolant s in PDRs ([Abel et al.|2 0051 [Holj 
lenbach fc Tielens||1997[ [Kaufman et al.||2006| ). We tried 
to tit the observ ed H2 line ratios to the FDK models by 
2006[ ) (see Figure [23|. The contribution 
H2 emission migHt impair the results. 



Kaufman et al. 
by shock excite^ 
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Fig. 11. — For the low resolution spectra the spatial evolution of several PAH and H2 lines over the whole slit for the ISOSS J20153+3453 
region i s pr esented. The IRAC 8 /xm map is shown and the 450 /xm SCUBA contours are overplotted. The slit overlay is explained in 
Section l4!6] 
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Fig. 12. Schematic representation of the ISOSS J20153+3453 
region. 

Nevertheless, the models show a trend toward higher gas 
densities in the outer layers of PDR (nn > 10^ cm~^) 
and UV radiation fluxes between 10^ and 10^ Gq. Sim- 
ilar results are found for ISOSS J04225+5150 East and 



GOlO.70-0.13 (Sections [5^ and [5^. Although these ob- 
servations can be interpreted as the result of a PDR, no 
indication of an associated H II region can be found in 
the 20 cm NVSS data. Unfortunately, there are no 6 cm 
survey products available that are sensitive enough for 
the detection of an H II region. 

One possible interpretation of the morphology of this 
region is shown in Figure [l2j We see the cationic PAHs 
and small grains photoelectrically heated by the IRAC 
sources at the western position. H2 emission and atomic 
fine structure lines originate from the warm layers of a 
PDR, but have also been formed via shocks. A clump 
of cold dust is located eastward in the foreground, ob- 
served in the submillimeter regime (Figure 10). The LL 



spectrum extracted at the eastern position is centered 
on the eastern MIPS 24 /im peak. Here a young, deeply 
embedded continuum source is driving the heating of the 
submillimeter clump. 



in the mid-infrared and submillimeter regime. The east- 
ern clump has a mass of M = 510 Mq and a temperature 
of Td = 17 K (Birk mann |[2QQ7| ). 

The IRAC images show a bright point source with ex- 
tended PAH emissions in the 8 /im channel, and dif- 
fuse H2 emissions at 4.5 /im. At 24 /im, the source is 
surrounded by warm dust. Slightly shifted northward 
from the MIPS 24 jam peak position the center of a sub- 
millimeter dust clump is detected by SCUBA (Figure 
14|. ISOSS J04225+5150 East shares several other fea- 
tures with the ISOSS J20153+3453 region: a class A 
PAH spectrum, [Si II] and [Fe II] lines. Hence, J-shocks 
and a PDR have been considered as the origin of these 
features as for ISOSS J20153+3453 West. However, in 
contrast to ISOSS J20153+3453, we detected the pres- 
ence of the [N e II] line. [Ne II] typically appears towards 
H II regions (Abel et al. 2005). The presence of the 
[Ne II] line can also be exp lained by shock excitation. 



Hollenbach fc McKee ( 1989 ) predicted the presence of a 
strong [N e IIJ line only for J-shocks with high shock veloc- 
ities. The [Ne I I] line ratios in Table Ig] agree with the pre- 
dictions from Hollenbach & McKee (1989). The H2 line 
ratios are similar to ISOSS J 20153+3453 West, therefore 
the results are similar for the PDR modeling (Figure 23). 
The line ratio fitting favor models with higher gas den- 
sity (nn > lO^cm"^) and UV radiation fiuxes between 
10^ and 10^ Go. For ISOSS J04225+5150, both PAH 
and H2 are present, with almost constant line strengths 
in the western direction. The source is associated with an 
elongated structure of higher dust density observed with 
SCUBA at 450 /im. The absorption feature at 10 jam 
originates from an outer ridge of dust. 




5.2. ISOSS J04225^5150 East 



ISOSS J04225+51 50 East 
8|jm 



5-10''AU 



Fig. 13. Overlay of the IRS slits for ISOSS J04225+5150 East. 
Notations are the same as those for Figure [9] 

ISOSS J04225+5150 East shows three compact sources 



Fig. 14. IRAC composite image of the ISOSS J04225+5150 
region: the 3.6 /im channel is colorized in blue, 4.5 /im in green 
and 8.0 /im in red. The contours of the SCUBA observations at 
450 /im are overplotted in yellow. The peak position of the MIPS 
source at 24 /im , which is used for the IRS pointing, is indicated 
by the red box point. 



5.3. GOlO.70-0.13 



[!h] 



The GOlO.70-0.13 refflon can be found in the Spitzer 
dark cloud catalog (Peretto fc Fuller|2009 ) . Two spectra 
were obtained for this region. The IKAC composite im- 
age shows dark filamentary structures against the bright 
PAH background. Around the eastern MIPS 24 /im peak 
position (GOlO.70-0.13 II), extended PAH emissions are 
present in the 8 /im band, as are warm dust emission 
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Fig. 16. Overlay of the IRS slits for the GOlO.70-0.13 region. 
Notations are the same as those for Figure [o] 

in the MIPS maps at 24 /im. A second western MIPS 
24 /im source is centered on a strong IRAC point source 
(GOlO.70-0.13 I). It is located in the tail of an extended 
submillimeter clump obtained from the SCUBA legacy 
program (see Figure 17). For GOlO.70-0.13 I, we iden- 
tified several absorption features, such as NHj, water 
ice, and CO2 ice, that are typical for deeply embedded 
sources. These are also observed for ISOSS J18364— 0221 
West as well as the sources in our sample associated 
with young, embedded sources (discussed in section 5.7 
). Even though no extended II2 emissions were found, 
the flux ratios for [Ne II] , [Si II] and [Fe II] (Tabl e [6|) for 



aiiuis 
ocated 



GQ1Q.7Q-Q.1 3 I+II agree with ratios for J-shocks 
et al.||2QlQD . The spectrum of GOIO. 70-0.13 II, 
to the west, is dominated by PAH bands. Bright UV 
radiation entails a class A PAH spectrum. 

Similar to ISOSS J04225+5150 East and 
ISOSS 20153+3453 this region shows indications 
for J-shocks and a PDR as well. The comparison of 



the H2 line ratios with the PDR in Figure 23 shows a 




Fig. 17. IRAC composite image of the GOIO. 70-0. 13 region. 
The 3.6 /im channel is colorized in blue, 4.5 /im in green and 8.0 /im 
in red. The peak positions of the MIPS point sources at 24 /im are 
indicated by the red box points. The IRS slits were centered on 
the position. 



correlation for a gas density of n > 10^ cm ^ and a UV 
flux of 1 - 2 X 10^ Go for GOlO.7-0.13 II. The line ratio 
fit for GOIO. 7— 0.13 I did not converge as for the prior 
spectrum; the eastern spectrum is likely just tracing 
the outer parts of the PDR, and the contribution of 
shocks to the H2 lines is stronger in this spectrum. We 
obtained the 6 cm and 20 cm survey products for this 
region f rom the Multi- Array Galactic Plane Im affln^ 
Survey ([Helfand et"aL| |2006) iWhite et al.| |2005D . No 



indications lor an associate^ 



White et al.||2005[ ). 

H II region were found. 
Our hypothesis is that GOIO. 70-0. 13 contains a PDR 
connected to a very young H II region, which is still 
optically thick at 6 cm and 20 cm. In such a case, the 
H II region would not be detected given the sensitivity 
of these radio observations. 

5.4. ISOSS J18364-0221 East 




Fig. 18. Overlay of the IRS slits for ISOSS J18364-0221 East. 
Notations are the same as those for Figure [o] 

Interferometric maps of 1.3 and 3.4 mm continuum 
emission show two cores. One core appears as a point 
source on the MIPS maps, surrounded by warm dust 
emissions at 24 /im. The IRAC observations do not show 
any point source at this position. The dust mass of both 
cores was estimated with 12 Mq and 18 M©, with a 
dust temperature of 22 K and 15 K. Narrowband NIR 
observations at 2.122 /im (H2 5'(1) 1-0) are related to 
outflow lobes, traced by C0(2-l) emi ssions around the 
MIPS source jHennemann et al.| |2009 ). 

The IRAC bands for 4.5 /im and 5.8 /im show some 
filaments, tracing H2 stemming from the outflow. The 
MIPS 24 /im source is surrounded by extended emissions 
in the 8 /im IRAC band. This band reveals a brighter 
eastward ridge. 

Although the SL orders only show a non-uniform ex- 
tended background spectrum, a strong PAH band ap- 
pears over the whole slit at 11.2 /im. The only detected 
lines are 5'(1) H2, [S I] , [Fe II] and [Si II] , which may 
be sho ck related. The [S I ] emission is not strong in 
PDRs ( [Kaufman et al.||2006|), but can be found, veloc ity 
independent, in J-shocks ( [Hollenbach fc McKee| |1989|) . 

Due to the lack of the SL order spectra, no further 
conclusion can be drawn. 

5.5. ISOSS 323053^5953 East 

The ISOSS J2 3053+5953 region has been studied by 
( |2007| ). It is located at a distance of 
region, two submillimeter cores with 



Birkmann et al. 
3.5kpc. In this 
Td = 19.5 K /17.3 K and a mass of ^ 200Mq each 
were observed. Both cores show signs of infall, as in- 
dicated by interferometric HCO+ measurements. IRAC 
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around the source position indicates that a predominant 
fraction of the gas is thermahy excited by the central 
source. The cooler gas traced by the 5'(1) line is ex- 
tended over more than 2.5^ Absorption bands of CO2 
ice and amorphous silicates indicate the presence of an 
outer ridge of cold material facing the observer. 

5.6. ISOSS J 18364- 0221 West - An isolated, icy source 



Fig. 19. Overlay of the IRS slits for ISOSS J23053+5953 East. 
Notations are the same as those for Figure [o] 

and NIR imaging show several embedded point sources, 
surrounded by diffuse PAH emission in the IRAC 4.5 /im 
images. Narrow-band observations (A = 2.122 /im) taken 
by the Omega 2000 prime focus camera at the Calar Alto 
3.5 m telescope show elongated H2 S{1) 1-0 shock exci- 
tation. 

The 24 /im MIPS observations show a bright source 
centered on the eastern submillimeter core, surrounded 
by extended warm dust emission and two fainter sources. 
Another faint source appears eastward, close to the peak 
position of the western submillimeter core. 

ISOSS J23053+5953 appears as a deeply embedded 
source. The extinction, calculated from the 9.7 jam sil- 
icate feature, is highest in our sample (see Table |5|. 
The bright 11.2 /im PAH feature is spatially extended 
over more than 1' with almost constant line strength, 
as shown in the spatial line fits (Figure 22). The spec- 
trum of ISOSS J23053+5953 East is dommated by neu- 
tral PAHs. No spectral indications for a PDR such as, 
e.g., strong cationic PAH bands were found in this re- 
gion. Furthermore, by fitting the H2 line ratios to the 
PDR models by Kaufman et al.| ([2006) no correlation 
between density n and FUV intensity G n m odeled for 
different line ratios was found (see Figure [231) • 

Therefore the observed forbidden lines of [Fe II] , 
[Si II] and [S I] do not originate from the hot layers of a 
PDR. These lines can be for med in the dense hot post- 
shocked gas of J-shocks ( jHollenbach fc McK ee 1989). 
Further more, the observed |!S I| line is never strong in 



PDRs (Kaufman et al. 2006), but can be observed as 



a velocity independe nt emission line in J-shocks ( |Hollen 



bachfcMcKee|l989). The presence of highe r H2 rotation 



lines is typical lor the post-shock gas phase ( Lahuis et al._ 
|£006). Therefore, a possible interpretation ol the spatial 
aistribution of the S{b) H2 line is that it represents the 
distribution of the cooling shocked gas, (the S{5) line was 
only detected sidewards of the MIPS 24 jam peak posi- 
tion). However, fluorescence pumping cannot be com- 
pletely ruled out as an excitation mechanism. The ob- 
served "green and fuzzy" features in the IRAC bands 
and narrow band imaging (A = 2.122 /im, H2, S{1)^ 
1-0) reveal the presence of excited H2. This supports 
the hypothesis of molecular hydrogen emissions from a 
post-shock relaxation region. The 6.22 /im PAH band is 
detected towards most of the positions where the S{5) 
H2 line was observed. Therefore, the excitation mecha- 
nism could be the same for this particular PAH as for the 
S{5) H2 lines. No molecular emission from CO, H2O, or 
OH, indicators for C-shocks, was f ound. Also the line H2 
ratios disagree with the results of [Kaufman fc: Neufeld 
(1996). The symmetric decrease of the S{3) line fluxes 




Fig. 20. Overlay of the IRS slits for ISOSS J18364-0221 West. 
Notations are the same as those for Figure [9] 
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Fig. 21. — Optical depths of the 15.2 /im ice feature for ISOSS 
J18364-0221 West. The solid line represents the low-resolution 
spectrum. The long dashed line represents the high-resolution 
spectrum. The two other lines represent laboratory absorbance 
data for an ice mixture of met hanol, water, and c arbon dioxide for 
temperatures of 70 and 80 K ( [White et al.|2009| >. 

The ISO SS J18364-Q221 re^on 2.2 kpc away has been 
studied by Birkmann et al. (2006). These observations 
show a large scale submillimeter clump with two separate 
cores in the SCUBA bands. For the western component a 
mass of M = 280tgo ^^^^ ^ temperature = 12.8 K 
were estimated. 

The IRAC observations show a bright compact object, 
some weak extended emission features and a faint point 
source eastward on the first airy ring. The bright point 
source on the northern edge of the airy ring of the central 
source must be treated as an artifact since it does not 
show the characteristic IRAC PSF. In the MIPS bands a 
single point source appears centered on the central bright 
IRAC source. The cross dispersion profile appears to 
be oblate and broadened, hence the target is consistent 
with the IRAC observations of a compact, but extended, 
source. The spectroscopy shows several ice absorption 
features and a broad silicate absorption at 9.7 /im. 

We used the IRAC and MIPS luminosities and the 
silicate extinction (Table [s]) with the SED fitting tool 
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Fig. 22.— Spatial line fits are the same as in Figure [TT] for the ISOSS J23053+5953 region. The IRAC image at 8 fim and the SCUBA 
contours at 450 /im are used. 
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(H2 0-0S(2) 12.3 /xm /H2 0-0S(1) 1 7 Aim ) 
(H2 0-0S(2) 12.3 Aim /H2 O-OS(O) 28.2 fim ) 
(H2 0-0S(1) 17 Aim /H2 O-OS(O) 28.2 A^m ) 




10^ 10^ 10-^ lO"^ 10^ 10^ 10^ 10^ lO"^ 10^ 10^ 10^ 10^ lO"^ 10^ 




10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 

n (cm-^) n (cm"^) 



Fig. 23. — Best fitting PDR models for certain H2 line ratios from [ Kaufman et al.| ( [2006| > are shown. The line fluxes are obtained from 
the high resolution orders. For visualization we used the PDR toolbox (http: //dustem. astro .umd. edu/pdrt/) . To show the contrast, 
ISOSS J23063+5963 is plotted as a region which does not host a PDR. 




5.9 6 6.1 6.2 6.3 6.4 6.5 7.2 7.4 7.6 7.8 8 8.2 8.4 8.6 8.8 15.5 16 16.5 17 17.5 18 18.5 



X [urn] X [urn] X [nm] 



Fig. 24. — Th e 6.2 jum (left) , 7.7/8.6 fim (center), and 16/18 fim (right) PAH bands. The global and local continuum has been removed, 
as proposed by |Peeters et al.| (2002). The 17.03 /xm feature is attributed to the H2 S{1) line. The PAH features for G19.27+0.07 were 
obtained from tlie background spectrum. 
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Fig. 25. — Overlay of the IRS slits for the IRDC sources. Notations are the same as those for Figure [9] 



from Robitaille et al. (2007) to obtain an estimate for 
the source parameters. FSirlSOSS J18364-0221 West, 
the modehng favors masses below 2 Mq and accretion 
rates between 10"^ - 3 x 10""^ Moyr"^ The mod- 
eled envelope and ambient masses are proportional to 
the accretion rates and range from 0.1 — lOM©. This is 
consistent with the scenario of a deeply embedded young 
low-mass protostar which could presumably evolve into 



a low-to-intermediate mass star. Birkmann et al. (20061 
had already classified this object as a low mass class 1 
object based on the NIK colors. 

The absorption profile for the 15.2 fim ice feature can- 
not be explained by pure CO2. We compared the profile 
of ISOSS J18364-0221 West with different a bsorption 
profiles fo r H2O+CH3OH+CO2 mixtures from [White et| 
al. (2009). The best fit was found for ice mixtures heated 
to 70-80 K (Figure 21). Therefore, the protostar is in 
fact the heating source of the ice, and not a background 
object. 

5.7. Young embedded IRDC sources 

G019.27+0.07 appears in the IRAC bands as an iso- 
lated point source at the northeast extension of an IRD- 
Cfilament. The background outside the filament and the 
shape of the filament change depending on the wave- 
length. This makes the spectral extraction for this source 
quite difficult. Therefore the source spectrum with a 
polynomial background fit and the nodded background 
spectrum is shown in Figure [7[ The 5.8/im and 8/im 
bands show the extinction of the IRDC against the galac- 
tic PAH background. The IRAC source is also de- 
tected as a point source in MIPS24/70. An extended 
submillimeter clump is observed at the position of the 




Fig. 26. IRAC composite image for G019. 27+0.07. The 3.6 /im 
channel is colorized in blue, 4.5 jim in green and 8.0 jim in red. 
The peak position of the MIPS points source at 24 fim is indicated 
by the red box points. The IRS slit is centered on the position. 
The SCUBA observation at 850 /im is shown as yellow contours. 

IRAC/MIPS source as denoted by the contours in Fig- 
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The clump can also be seen a t 1.2mm, and the 
core mass was calculated as IISMq (jRathborne et al. 



2006 ). H2O and thermal CH3OH maser emi ssions were 
observed wi th the Green-Bank- Telescope by [Chambers 
et al. (2009). The later maser emission indicates a more 
evolved evolutionary stage. The background spectrum 
of G019.27-h0.07 reveals a A spectrum. These PAH fea- 
tures are not detected toward the source spectrum, hence 
the UV excitation source can be external. Molecular hy- 
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drogen emissions are found in the source and background 
spectra. This is consistent with a sUghtly extended blob 
of H2 observed in the IRAC band at 4.5 /im. 

G034.43+0.24 is an isolated source found in MIPS and 
IRAC images. The peak of the submillimeter clump at 
850 jam is centered on the 24 jam source position. The 
extended lobe in the 3.6 jam map of molecular hydro- 
gen lines could be excited by an outflow. The same 
absorption feature as that in ISOSS J18364-0221 West 
is observed. This supports the hypothesis of a young, 
deeply-embedded object with outflows. For the remain- 
ing targets embedded in IRDCs, no extended hydrogen 
features were resolved on the IRAC maps. However, 
molecular hydrogen emission was detected in G028.04- 
0.46 and G053.25+0.04 I+II. 

5.8. Comparison with Orion IRc2 

Many of the general features observed in the ISO spec- 
tra of our sample are also observed in the spectrum of 
Orion IRc2 (see Tabled. 

The 9.7 /im silicate absorption feature can be found in 
almost every source. As in IRc2, we find ice features of 
water, methane and CO2. NH4 in the gas phase appears 
in absorption. Several molecular features from IRc2, e.g., 
SO2, C2H2, and CII4 (in the gas phase), are not present 
in our sample at the IRS resolution and sensitivity. Most 
sources show a wide variety of PAH bands. Since the 
times of the ISO observations of IRc2, there has been re- 
mar kable_2rogress in the understanding of PAH composi- 
tion ( Tielens|2008 ). We constrained the physical proper- 
ties ot the tar- U V for several sources from the PAH bands 
at 6.8, 7.7, and 8.6 /im. Due to the lower spectral resolu- 
tion of Spitzer IRS compared to 750 /SWS, features like 
H2O in the gas phase can not be observed in our sample. 
Since the spectra of IRS are taken longward of 5.5 /im, 
only the 0-0 5'(0)-5'(7) lines can be observed. The Q and 
O modes and the higher orders for S remain inaccessible. 
It is difficult to distinguish between the different excita- 
tion mechanisms for molecular hyd rogen, as was done for 
the Orion peaks ( Rosenthal et al.|[2000) , because we are 
missing these orders at shorter wavelengths. However, 
the spatial distributions of the 5'(5) and 5'(3) lines en- 
able us to discriminate between thermal excitation and 
shock excitation. For some of the sources in our sample 
the [Ne II] line is observed in our sample as well, and 
can be correlated to a PDR and/or J-shock excitation. 
In contrast to the Orion star forming region we did not 
find any indications for C-shocks. All spectra in our sam- 
ple show only a fraction of Orions IRc2 characteristics. 
Either they show deeply embedded young sources with 
strong silicate and ice absorption bands, or more evolved 
sources with indications of PDRs and strong PAH bands 
below 9 /im. 

6. SUMMARY 

We used the IRS on board of the Spitzer Space Tele- 
scope to observe a sample of 14 intermediate-mass YSOs 



in different evolutionary stages. We present low- and- 
high resolution mid-infrared spectra of these sources. 
The IRS spectra demonstrated how difficult it is to dis- 
entangle the complexity of such starforming regions, be- 
cause the IRS beam covers multiple physical components. 
However, with the help of previous mult i- wavelength 
studies we could shed light on several physical compo- 
nents and properties. In contrast to previous spectro- 
scopic ISO observations, we derived information about 
the spatial extent of lines using the IRS slit spectroscopy. 
The spatial distribution of selected PAH and H2 lines in 
comparison with previous imaging observation can shed 
light not only on the YSOs but also on their vicinity. 
The infrared spectroscopy shows two different kinds of 
sources: 

• Young isolated sources, which are still embedded 
in an envelope of cold dust and ice. These sources 
tend to still have low to intermediate masses with 
high accretion rates. Only neutral background 
PAHs are found for these sources. Several different 
absorption features, such as NHj (6.75 /im), wa- 
ter, methane, and CO2 ice, can be found for these 
sources. CO2 ice can also be found for the more 
evolved sources too. For ISOSS J23053+5953 East 
the same indications for J-shocks are found as for 
the more evolved source ISOSS J20153+3453 West 
(see below). 

• The second group of sources is more evolved. We 
found H2 pure rotational lines, atomic fine struc- 



ture lines such as [Fe II 
bands. The presence of 



[Si II] and several PAH 
Si II] , [Fe II] and H2 lines 



can be explained by J-shocks. For some sources 
[Ne II] indicates high- velocity shocks. The line ra- 
tios for [Ne II] and H2, [Si II] , [Fe II] agree with the 
predictions for J-shocks. Shock excitation might 
also explain the spatial distribution of certain H2 
lines for ISOSS J20153+3453 West. There are no 
indications for C-shocks in our sample. For some 
sources a class A PAH spectrum was detected with 
strong PAH bands below 9 /im. All these sources 
have indications for a PDR. Although the contri- 
bution of shock excitations t o the H2 line eniission 
is not known, PDR models ( Kaufman et al.|[2006 ) 



favor higher gas densities (> 10^ cm and UV 
field strengths between 10^ and 10^ Go- 
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